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ABSTRACT 


Whittaker's  solution  of  the  wave  equation  is  written  in 
a  form  from  which  all  solutions  used  in  problems  of  propagation 
of  a  disturbance  may  be  derived.  An  elementary  theory  of  so- 
called  "refraction"  waves  clarifying  the  physical  conditions  of 
their  generation  is  discussed. 
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PART  I 

On  the  General  Solution  of  the  Wave  Equation 

1 .  Introduction 

Although  the  propagation  of  a  disturbance  produced  during 
an  earthquake  is  a  special  case  of  a  general  problem  of  wave 
motion  in  a  continuum,  there  is  a  specific  difficulty  arising  from 
the  fact  that  the  nature  of  the  source  of  this  disturbance  remains 
hidden  in  the  Earth's  interior. 

In  many  investigations  following  Lamb's  (1904)  well  known 
paper  attempts  were  made  to  find  the  solutions  representing  the 
disturbance  propagating  from  a  given  source  in  each  particular 
problem.  The  methods  used  in  these  investigations  diverge  on 
several  points  but  the  solutions  of  most  problems  are  built  up 
from  a  set  of  particular  expressions  which  are  generalized  until 
all  the  conditions  of  a  problem  are  satisfied.  It  is  obvious  that  a 
more  logical  and  natural  way  would  be  to  make  use  of  the  general 
form  of  solutions  of  equations  of  the  wave  theory  and  by  specifi¬ 
cation  to  obtain  those  expressions  which  will  satisfy  all  the  con¬ 
ditions  of  a  particular  problem. 

The  propagation  of  a  disturbance  in  a  medium  is  determined 
by  one  or  more  functions  (components  of  a  displacement  in  an 
elastic  medium  or  components  of  the  electric  and  magnetic  force) 


3. 


which  must  be  solutions  of  one  partial  differential  equation  or 
of  a,  system  of  such  equations.  A  complete  integral  of  these 
equations  is  used  in  b  oundary  problems  and  the  number  of  arbitrary 
constants  should  be  such  that  they  all  can  be  determined  by  means 
of  boundary  conditions.  On  the  other  ,hand  it  may  be  important  to 
begin  an  investigation  with  the  general  integral  of  the  given  system 
of  partial  differential  equations  and  to  make  it  clear  which  physical 
conditions  in  z  particular  problem  require  the  choice  of  some 
special  functions.  Such  an  approach  to  the  solutions  of  fundamental 
problems  of  Seismology  is  the  principal  purpose  of  this  paper. 
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2.  The  general  solutions  of  the  wave  equation.  In  order  to 
determine  the  propagation  of  ,a  disturbance  in  an  elastic  medium 
we  have  to  .consider  three  equations  of  the  form 


(i) 


which  must  be  satisfied  by  the  components  (u,  v,  w)  of  a  displace¬ 


ment  S  ,  if ^  is  the  constant  density  of  the  medium,  $=divS  and 
X  and^tare  the  elastic  constants.  The  three  unknown  functions  may 
be  expressed  in  a  simple  form  in  terms  of  four  others  %  Y'.Yi’K 
i.  e.  in  terms  of  a  scalar-  and  vector -potential.  Then,  on  putting 


yicccicf  +  Cial 3  ,  fiCYu'fiL,n)  5 


(2) 


^ Z 


and 


i  X  + 


The  equations  of  motion  (1)  are  satisfied  by  the  expressions 
and 'f/- are  solutions  of  two  wave  equations 


-  d* 


(2),  if  f 


(4) 


Thus,  (4)  represent  a  system  of  sufficient  conditions  for  the 
equations  (1).  It  ?eems  that  it  has  not  yet  been  proved  that  (4)  are 
also  necessary  conditions.  We  also  add  a  second  remark.  Because 
of  the  transformation  (2)  a  completely  arbitrary  element  is  intro¬ 
duced  as  a  fourth  function.  Therefore  a  further  condition  is  usually 
imposed,  for  example  ,  In  a  simpler  case,  where  there  is  a 
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plane  boundary  in  a  medium  (or  boundaries)  and  the  propagation 

can  be  described  in  terms  of  cylindrical  coordin  ates  (r,  z),  only 

two  components  of  displacement  are  considered.  These  components 

are  expressed  in  terms  of  two  functions  and^,  i.  e.  the  equations 

similar  to  (2)  introduce  only  two  unknown  functions. 

The  boundary  ^.nd  initial  conditions  must  also  £e  given  in 

each  particular  problem  of  propagation  of  a  disturbance,  but  we 

/  genera  (?  So ~t'&n 

will  now  write  different  forms  of  thetwave  equation.  On  making 
use  of  a  polynomial  linear  in  coordinates  and  time  and  of  two 
parameters  U,  and  the  general  solution  of  (4)  was  giver ,  by 
Whittaker  (1947)  in  the  form 

siniL$inir-hZco$  u+oct,  u,v)afa  dir  (5) 


f/j  fo 

-n  drr 


z  Si'n  U  COS  V“-h 


1 


and  by  Bateman  (1904),  with  a  single  parameter  0  , 

<f=  I  ■F(xces4+vsuifi+iZ,  y+t2*n6  +  cctc«se,  d& 


(6) 


*// 


where^ and^  are  arbitrary  functions.  It  is  easy  to  verify  that  these 
functions  satisfy  (4).  It  is  also  well  known  that  Lamb,  Sommerfeld, 
Weyl  and  other  investigators  £>egin  with  a  special  form  of  the  function 
namely,  the  exponential  function. 

Whittaker's  solution  as  a  general  integral  is  expressed  in 
terms  of  one  arbitrary  function  and  does  not  depend  on  some 
arbitrary  constants,  the  integrations  with  respect  to  two  parameters, 
u  and  v,  being  performed.  As  to  the  existence  of  Bateman's  solution, 
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it  is  known  that  there  may  be  different  forms  of  the  general 
solution  of  a  partial  differential  equation  which  cannot  be 
transformed  into  one  another.  In  the  following  discussion  use 
of  (5)  only  will  be  made. 

It  is  obvious  at  once  that  the  integrand  in  (5)  can  be 
written  also  in  the  form 

f  +  oct  9  (7) 

since  the  first  equation  (4)  is  satisfied  by  (7)  provided 

ajL-hS*'+ct=  /  (8) 

Thus,  for  example,  we  can  put 

- CSih.fi  u  CiTShir  9  S--Sihk  USffih  irt  c  -cosh  1C  (9) 


and  write  instead  of  (5) 

OO  JO 


—  OO  —  00 

The  number  of  parameters  in  Whittaker's  solution  may  be 
increased,  for  example,  on  writing 


««  OJOfl 


V-...+  V.t  +  w,  VVlW)du^ATliW(l\) 


-09  -OO  —00 


and  it  seems  to  be  natural  to  take  the  argument  representing  the 
most  general  form  of  a  linear  function  of  four  independent  variables, 
An  arbitrary  function  ^  of 

/v 

jl  =  iffx-*;  +  rt  (y  -$)+■  v3  (z-z)+  (* ' t)  (12> 
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A  A  0 


depending  also  on  parameters  Yty  will  satisfy  the 

equation  provided 


wave 


(13) 


Then,  the  expression 

OO  oO 


f'J-J  {( 

-<*>  -oo 

is  the  corresponding  form  of  the  general  solution  of  (4),  if  all 
special  mathematical  conditions,  like  the  convergence  of  integrals 
etc.  ,  are  satisfied.  Moreover,  because  of  (13)  we  can  make  use  of 

two  arbitrary  functions^  and  ^depending  on  two  different  expressi 

\ 

for  the  function  (12),  namely, 

f-z~V,  (X-X)  *  Vx  (y-p  -  vz(z-  Z )  ^Cb' 


sions 


(12') 


Thus,  we  can  write 

oo 


? *  P’J i Ch ;  )<**,-.<&  + 


(14') 


-OO 

oa 


-t- 


P-\  {*-(h  ’  •  ••  ^ 

-  oo 


As  a  special  case  we  can  consider  a  function 

LA  (n  e,*) 


<f=T.eL> 


(is) 


where  P  is  an  arbitrary  function  of  parameters  .  As  to  the 


limits  of  the  integrals  in  (14)  there  can  be  certain  reasons  to  assume 
them  to  be  equal-^a^  or  to  jhave  some  others.  In  general,  if  we  want 
to  apply  the  general  solution  in  the  form  (14)  to  pome  special  problem, 
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we  have  to  attempt  to  adjust  all  arbitrary  elements  involved  to 
the  conditions  of  the  problem. 

As  to  the  complete  integral  of  the  wave  equation  it  is  easy 
to  see  that  it  should  depend  on  four  arbitrary  constants.  In  the 
general  case  of  one  partial  differential  equation  of  the  second 
order  with  four  independent  variables,  if  a  function  . tfX.^  f 

Cn...CK)  is  a  complete  integral  it  depends  on  K  constants,  K  being 
the  number  of  the  first  and  second  derivatives  of  2  with  respect  to 
X  ,  i.  e.  14.  Since  there  should  be  certain  relations  between 
these  constants  satisfied  if  a  differential  equation  has  a  simpler 
form,  the  number  of  independent  constants  will  decrease.  It  is 
obvious  that  there  is  a  certain  function  z(pC,}..  (^depending  on 

4  constants  which  will  yield  a  wave  equation. 
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t 


3.  Some  special  forms  of  the  solution  (14).  We  are  able  now  to 
show  in  which  way  the  known  forms  of  solutions  follow  from  the 
expression  (14)  by  a  specialization  of  conditions. 

1.  -Point  source  in  a  homogeneous  mechurm  If  this  source  is  at  a 
point  S(0  ,  0  ,  h)  we  have  to  put Z~rJ.  If  the  source  begins  to 
emit  the  disturbance  at  t^O  we  put  also  t"0.  It  is  easy  to  see  that 
we  have  to  make  use  of  the  function  (15)  when  S  is  a  source  of  an 
harmonic  oscillation.  This  function  now  takes  the  form 

y„  )  e  '  (i6) 

provided  the  condition  (13)  is  satisfied.  For  harmonic  waves  the 
angular  frequency  Despite  the  fact  that CJ  is  an  essentially 

positive  quantity,  we  can  make  use  of  the  limits  CJ-+00 ,  since  for 
P  even  in  Vy  and  real  we  can  use  the  function  Cos  Y^t.  and  for  P  odd 
in  Vy  the  term  P$ irvV^'t .  The  integral  from-**  to-'-co  equals,  then, 
tv/ice  the  integral  from  0  to  oo. 

The  general  solution^? now  takes  the  more  simple  form 

f=j  j’j'  '?■(%,  ,  Xy  )  +  (1?) 

-  (to  — OO  -  vo 

if  the  source  is  supposed  to  emit  all  frequencies  from  0  to  vO 
(ox  -YsCZ-^J  in  exponent). 

On  the  other  hand  we  can  assume  that  Y^  ~  - .  Then,  the  time 

q:  pt 

factor  takes  the  form  C  '  ,  (if  p  is  real  and  positive,  this  factor 

corresponds  to  a  disturbance  continuously  decreasing  or  increasing 
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at  each  point.  The  solution  corresponding  to  the  last  case  was 
considered  by  Cagniard  (1939). 

The  condition  of  spherical  symmetry  holds  for  a  point 
source.  Having  in  mind  further  applications  we  consider  symmetry 
with  respect  to  the  z-axis.  Let 

cc-  z,  ces  ic  7  vt  -  jt'  cos  r  ?  K  [d,  o°  )9 

y  *  r  u*  }  *  sear  ,  zr  (~/ry  n ) . 


(18) 


Then,  (17)  takes  the  form 


r 


oo 

=Je 

-  o o 


v;)«rWr 


(19) 


Since,  by  (13)  and  (18), 


(20) 


and  does  not  depend  on? 

,v>0 


we  can  write 

=  /  ecy,,tcM  / e  ^ (Z  r,  v^e 


iTHCob  (i  -  h) 

u>  =  f  e  cW¥  /  «  an  j  dz  (21) 

-  OO  0  -n 

If  the  medium  is  an  isotropic  one,  there  is  no  reason  to  have 


in  (21)  a  factor^  depending  on  X.  and  u.  Therefore,  we  can  use  the 
simpler  expression 

eo  rOO  r 'T 

(  cv^t  ,  /  LV,  (Z-k)zz  i  .  /  \  /  /  / 

Cf  =  J  e  dv¥ /e  *  yr  (22) 

-  oo  0  ~7i 

The  last  integral  represents  the  Bessel  function  (See,  for 


example,  Heine  (187 8)j 


K 


J0(*z)=zle  d 


71 

Ie 


2 


(23) 


-/* 


whent£is  an  arbitrary  term.  Thus,  the  factor  1/2H*  being  included 
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into  the  arbitrary  function 
c?0  /-  G0 


f  ,yt,  fZr  T 

<f~Je  *  J  <Ak. 

_  ix>  0  _ 


(24) 


If  we  make  the  assumption  that  ^ Q^. ;  (Y^j, 


<*>  ■  ^  r*00 

Y9f'Z(v*')e‘r*  ^  frfo  Jo(lCT)e 


(25) 


/) 

Now,  it  is  easily  seen,  that  in  order  to  obtain  the  well  known 


expression 


■  #>  -y/2-Af-  , 

&  J0Lfir )  *<**■ 


(26) 


where  f\.  ,  corresponding  to  spherical  waves  emited  by  one 


source,  we  have  to  put 


y=+^i  , 


(27) 


the  upper  sign  has  to  be  taken  if  z  -  h  is  positive.  Finally 


e  cUo 


(28) 


if  ^Lf-CO  .  The  function^^^ojis  the  last  one  not  yet  fixed.  The  most 
natural  assumption  is  that  ^  is  connected  with  the  properties  of  the 
source.  The  restrictions  concerning  the  function  ?  in  (21)  is  now 
made  clear  by  the  expression  (26);  these  restrictions  are  necessary 
to  obtain  the  potential  Cp  cor  re  sponding  to  a  spherical  wave. 

In  all  expressions  the  Riemann  integrals  were  used,  but  it  is 


obvious  that  Stieltjes  integrals  could  also  have  been  used.  The 
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integration  with  respect  to  the  parameters  in  (14)  is  not 

subject  to  any  special  restrictions  and  it  seems  that  there  are  no 

( 

objections  to  the  use  of  the  integration  over  the  cj  or  complex 
plane.  This  is  done  by  Lapwood  (1949)  under  the  assumption  that  the 


initial  disturbance  has  the  form 

/ 


Jl 

where  si  denotes  a  strip  of  the  ^complex  plane  below  the  real  axis. 


t<  0 
tr>0 


Thus,  Lapwood' s  solution  follows  from  (14)  or  (25)  when 


/SCO 


and  the  integration  is  changed  into  that  over  the  strip  jCX  . 

It  is  well  known  that  if  a  source  produces  a  disturbance  in  a 
bounded  homogeneous  medium  the  direct  wave  is  given  by  the  ex¬ 
pression  (28)  and  the  total  disturbance  at  a  point  P  (r,  z),  z  ^  0 , 
at  a  time  t,  may  be  represented  by 

(29) 

f 

The  term  ,  if  the  boundary  is  a  plane  z  =•  0,  may  be  taken 
in  the  form 

S®  ^  m  oo  -y 

c/^J?t(K)J(n)e‘  }  dK  ,  z>o.m 

•  oe  0 

Niessen  (1933)  has  shown  that  a  function  of  this  kind  (time 
factor  b  eing  omitted)  is  a  sum  of  two  terms,.  The  first  represents 
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the  reflected  wave  and  the  second,  written  in  the  form  of  an 
integral  over  the  plane  z=0,  suggests  the  action  of  points 
along  the  boundary  in  the  sense  of  the  Huygens  principle. 

It  is  obvious  that  a  boundary  surface  other  than  a  plane 
would  require  a  more  complex  form  of  function  . 

The  solution  (25)  corresponds  to  the  assumption  that  an 
arbitrary  initial  disturbance  can  be  always  built  up  from  a 
harmonic  solution,  i.  e.  is  the  so-called  steady-state  solution 
used  by  most  investigators.  Now,  equation  (25)  follows  from  the 
general  integral  of  the  wave  equation  and  by  the  Fourier  formulae 
we  have  :  if 

oO 

icet  iff  Os  t 

-oa  -Do  -oo 


On  comparing  this  equation  with  (25)  we  obtain  the  function 
on  putting  Y ^  -  Co  and 


Thus,  the  solution  (14)  was  reduced  to  this  well  known  form  because 
of  the  special  choice  of  the  function  f  ,  given  by  equation  (16).  The 
single  remaining  arbitrary  factor,  the  function  'P  (V,  can  be  used  in 
order  to  satisfy  the  boundary  conditions,  as  is  well  known. 
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2.  Line-source.  If  there  are  several  sources  acting  on 
the  z-axis  at  the  points  z=hp  h^,  ...»  hg  we  should  use  the  sum 
of  functions  of  the  type  (25)  provided  all  other  conditions  are 
similar  to  those  mentioned  above.  Coulomb  (1931)  considered 
continuous  distribution  of  sources  on  the  z-axis.  If  we  assume 

A  s\ 


that  in  (14)  X  ~^j~0  and  write 


(33) 


this  expression  would  correspond  to  the  case  where  there  are 

sources  at  all  points  between  z  *  and  z  It  is  obvious  that 

Coulomb's  generalization  . 

-/7c  71 

o  A 

where  R^*  and  Z,  our  notation,  is  a  special 

case  of  (33). 


3.  Sources  on  a  surface  or  in  a  volume  placed  at  a  finite 
number  of  points  ZsViU  produce  a  potential 


For  a  continuous  distribution  we  have 


(35) 


(36) 


where  D  is  a  two  or  three  dimensional  region.  In  the  first  case  the 
last  integral  must  be  transformed  into  a  double  one,  X  -  *  ) 
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being  the  surface  along  which  the  sources  lie. 


On  assuming  Z-0  and  on  taking  again  the  integrand  in  (14) 
in  the  form  (16)  we  obtain  the  expressions  for  the  scalar  potential 
Cp  and  the  vector  potential ^Jwhich  were  used  by  Schermann 
(1946).  He  puts 

i[y,  (*-*  i+iQf  -p '  ^ 


Sv 


TJ 


-00 


'Pe 


and  similar  expressions  for  y^-  in  order  to  obtain  the  solution  of  the 
problem  of  oscillations  in  a  semi-infinite  solid  when  the  displace¬ 
ments  or  the  external  forces  at  the  boundary  are  known.  Hallen 
(1936)  used  a  similar  form  in  a  two-dimensional  problem.  Our  form 
of  the  general  solution  (14)  of  the  wave  equation  make  it  clear  at 
once  that  (37)  may  be  in  terpreted  as  a  potential  produced  by  sources  at 
the  boundary  plane  z  ro.  On  applying  his  method  Schermann  meets 
the  difficulty  mentioned  in  paragraph  2  concerning  the  representation 
of  three  variables  u,  v,  w  in  terms  of  four  functions  such  as  ^  ^ 

There  are  also  four  arbitrary  functions:  P  in  the  expression  (37) 
and  in  similar  ones  written  for  ft  If,  for  example,  the  displace¬ 
ments  a*  the  boundary  plane  are  given  functions  of  x,  y  and  t, 

P  and  have  to  be  expressed  in  terms  of  ~U0  ,  V0)w0  .  Since  there  are 
only  three  conditions  which  we  obtain  on  substituting  ir^wm)  and 
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similar  expressions  for  into  (2)2,  (X  can  have  an  infinite 
number  of  different  values.  Schermann  takes  "the  most  simple 
ones".  A  physical  argument  would  be  necessary  in  this  case. 


4.  A  source  in  one  of  two  semi-infinite  media.  We  will 
now  discuss  several  other  points  concerning  this  problem  which 
we  began  to  consider  in  1.  /.  We  obtained  (25)  from  (14)  under 
several  assumptions  and  one  of  them  was  ~t  -  0  .  If  we  suppose  now 


that  ~t  differs  from  zero,  we  can  write 

Y,X  +  (13') 


The  expression 

<x> 

_  oo 


(38) 


represents  a  special  form  of  a  general  solution  of  the  wave  equation 
with  a  reduced  number  of  parameters,  which  should  correspond  to 
the  existence  of  a  single  source  located  at  the  point  (0,3  0,  h).  Moreover, 
if  we  assume  that  4  V^  -  p~  is  a  real  and  positive  number,  we  obtain 


Since  the  function  f  itself  satisfies  the  wave  equation  the 
double  integral  with  respect  to  \  and  is  also  a  solution  and  we  can 
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also  have 


f 


/» 

t. 


5  ^  t-tj 


(40) 


-bO 


where,  because  of  the  fact  that  p  and  t  are  here  assumed  to  be 

constant,  the  function  will  satisfy  the  wave  equation. 

From  the  formal  viewpoint  there  are  no  objections  to  summation 

A 

for  a  certain  interval  of  the  variable  parameter  t.  For  example, 
the  function 


•  A 

is  obviously  a  solution  of  the  wave  equation  if  4  and  ^are  two 
constant  values  of  the  parameter  t.  .These  can  be  0  or  -oo  and 
*MX>  ,  respectively,  and  it  is  a  natural  question  to  inquire  into 
the  meaning  of  a  summation  with  respect  to  the  parameter  t. 

It  is  obvious,  that  if  there  is  a  set  of  sources  placed  at  different 
points  which  begin  to  produce  disturbances  at  different  times,  then, 
in  the  expression  of  the  potential  the  corresponding  values  have  to 
be  ascribed  to  the  parameter  t.  This  may  be  generalised  for  a 
continuous  distribution  of  sources,  but  it  is  hard  to  see  any  signifi¬ 
cance  in  such  a  summation  in  the  case  of  a  single  source. 

If  we  denote  t  by  -u  and  assume  that  the  condition  of  axial 
symmetry  is  already  taken  into  account,  the  integrand  in  (41)  may 
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be  written  in  the  form  e/’w  ¥  (t,2,t+2c) . 

This  is  the  function  considered  by  Cagniard  (1939),  who 

A 

assumes  that  the  integration  holds  between  the  limits  2 — 00  } 


i.  e.  that 


^,2,^  Il)gLcs 


(42) 


also  represents  a  solution  of  the  wave  equation.  Whether  or  not  we 
can  use  variable  limits  for  the  parameters  in  an  expression  like 
(14)  deserves  investigation. 

There  is  no  doubt,  however,  that  the  formulae  (29)  and  (30) 

may  be  used  in  order  to  represent  the  scalar  potential  in  a  medium 

containing  a  source  and  that  we  can  assume  for  the  second  medium 

the  potential %in  the  form 

vo  oO 

J '^(*)Jt>CKZ)elV}Z<6c  ,  z<<?.(43) 


-  00 


0 


To  obtain  the  functions  and  as  well  as  the  components  of 
the  vector  potential  from  the  boundary  conditions  the  full  expressions 

(30)  and  (43)  are  not  needed  and  the  "reduced"  form  of  them 

+  ,  £=e‘y^f>v,*,z;  (44) 

can  be  used.  Written  in  the  form  of  integrals  which  can  easily  be 
seen  in  the  equations  given  above,  the  functions  X ^( >y , 2T;2T),  scud 
X2(Vv,X)2)  and  YjO Y2(Vy)t,Z)  corresponding  to  ,  the  "^-functions 
should  represent  all  kinds  of  waves  (the  direct  excluded)  which  may 
be  generated  in  both  media.  The  evaluation  of  these  integrals  presents 
well  known  difficulties.  The  divergent  results  obtained  by  Sommerfeld 
and  Weyl  in  the  analogous  problem  of 
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electromagnetic  theory  may  be  quoted  as  an  example.  Muskat 
(1933)  seems  to  be  the  first  who  gave  the  approximate  solutions 
in  such  a  form  that  it  was  possible  to  see  the  existence  of 
different  wave  fronts  travelling  in  both  media.  He  considered 
liquid  as  well  as  solid  media.  Another  method  was  followed  by 
Cagniard.  We  recall  thatil^is  p  in  Cagniard’ s  notation.  Now, 
it  may  be  assumed  that 


X(p,  X,x)= />J e.~nACs^'c^oLl’  (45) 

o 

y  (p,r,z)c/>J  (46) 

0 


From  these  equations,  four  other  functions  Aj,  B  j  and  A£,  for 
two  media  may  be  determined.  They  are  Laplace  transforms  of  the 
functions  X^/p  and  Y^/p  or  the  solutions  of  the  Carson's  equations 
(45-46),  in  which  the  variable  5  has  the  meaning  as  follows. 

It  was  shown  by  Cagniard  that,  if  we  put  the  parameter  k=pq, 
q  being  a  new  variable,  the  parameter  s  will  be  connected  with 
other  quantities  by  the  equation 

|  ^  y  j 

^  k  ^  +  z  =  s 

where  0^  and  ^  are  velocities  in  the  first  medium  and  C (or  more 
exactly  /T-6*'  )  is  that  used  in  (23).  The  investigation  of  solutions  of 
(45)  and  (46)  as  well  as  solution  obtained  by  other  methods  proves  the 
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existence  of  several  wave  types  in  both  media.  These  results 
will  be  discussed  in  the  second  part. 
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PART  II 

On  Waves  Generated  At  An  Interface 

1.  Introduction. 

In  order  to  explain  arrivals  observed  in  seismic  refraction 
shooting  "refraction"  paths  and  waves  associated  with  these  paths 
(refractions)  are  assumed.  These  waves  propagate  in  layered  media 
along  with  the  familiar  direct,  reflected  and  refracted  waves. 

The  experiments  of  O.  Schmidt  (1939)  provide  a  clear  demonstration 
of  four  types  of  waves. 

Refractions  have  the  following  characteristic  properties 
(figures  1  and  2):  (1)  They  are  represented  by  straight  lines  (W1^) 

joining  points  of  the  interface  reached  by  the  refracted  waves  (W2) 
and  points  T  at  which  they  are  tangent  to  the  reflected  wave  (  H^); 

(2).  If  the  source  is  placed  at  the  interface  as  in  Figure  4,  refractions 
are  tangent  to  the  direct  wave  (W^);  (3)  They  are  inclined  to  the  inter¬ 
face  at  an  angle  equal  to  that  of  total  reflection;  (4)  the  waves  (  W^) 
propagated  in  a  medium  of  lower  velocity  precede  the  reflected 
wave  and,  under  certain  conditions,  the  direct  wave. 

An  explanation  of  refractions  is  suggested  by  investigations  of 
a  steady  state  solutions  of  the  propagation  of  disturbance  in  layered 
elastic  media.  Jeffreys  (1926,  1931)  seems  to  be  the  first  to  obtain 


a  term  which  could  be  interpreted  as  a  wave,  which  "appears  to  have 
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travelled  along  the  interface  with  the  velocity  of  sound  in  the 
lower  medium".  Muskat  (1933)  giving  the  approximate  solution 
by  methods  of  Lamb  (1904)  and  Sommerfeld  (1909)  for  the  case 
of  longitudinal  waves  in  two  semi -infinite  layers  and  in  more 
general  case  of  two  superposed  elastic  media,  also  obtained  terms 


of  the  form 


A  -cwCz/vz  +  €  ft  +  z-) /v,  ) 


=  — ,  e 


(i) 


where:  A  is  the  factor  depending  on  elastic  constants  and  on  the 
frequency  a)  of  a  harmonic  disturbance,  r  and  z  the  cylindrical 
coordinates,  h  the  distance  of  the  source  from  the  interface  and 
vi<  V2  the  velocities  of  propagation  of  the  disturbance  in  the  upper 
and  lower  medium  respectively,  /=  O-  *!*/%  )*. 

Sobolev  (1932)  obtained  wave  fronts  similar  to  those  re¬ 
presented  by  the  lower  part  of  Figure  2,  if  we  assume,  namely,  that 
only  the  halfspace  2  is  an  elastic  medium  and  O  is  the  source  of 
compressional  and  distortional  waves. 

Since  that  time  the  "refractions"  were  treated  in  many  in¬ 
vestigations  and  were  given  different  names,  namely,  the  "fourth"  ray  and 
the  "travelling  reflection"  (Schmidt,  1939),  the  "conical"  waves 
(’Cagniard,  1939),  Minthrop-  or  head-waves  (Kopfwellen,  in  a  sense 
different  from  that  applied  to/4lhell)  or  side -waves  (Flankenwellen) 


by  seismologists  in  Germany  and  USSR. 
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In  the  Figure  2  are  represented  different  wave  fronts, 
which  were  given  in  this,  most  complete  form  by  Cagniard 
(1939,  p.  122).  They  correspond  to  the  case  of  two  semi-infinite 
elastic  media  and  one  point  source  of  compressional  waves,  S. 
These  waves  are  found  under  condition 

fi,  <  «,  *  h  <  **  <2) 

on  solving  the  integral  equations  (45) -(46)  Part  I.  They  are: 

Wj  -  the  direct,  -  the  reflected  compressional  wave,  -  the 

reflected  distortional,  W2  -  the  refracted  compressional  and  ^ 
the  refracted  distortional  waves.  All  other,  "conical"  waves 

1  ff  ,  ft 

(  ^  ,,,,  )  are  of  the  "refractions"  type. 

Merten  of  the  Shell  organization  seems  to  be  the  first  who 
gave  (1927)  the  correct  idea  about  the  generation  of  "refractions", 
using  the  "wave-front  diagram",  according  to  Thornburgh  (1930) 
(see  also  Ansel  (1930)  and  Dix  (1939)). 
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2.  The  elementary  explanation  of  waves  of  the  fourth  kind. 

The  characteristic  physical  conditions  under  which  the 
refractions  are  generated  may  be  obtained  in  a  most  simple  case 
of  two  semi-infinite  homogeneous  elastic  liquid  media  (Fig.  1). 

We  assume  that  a  spherical  compressional  wave  is  emitted  by  a 
source  S  located  at  a  distance  h  from  the  interface  zsO.  When 
this  wave  reaches  the  interface  at  the  base  of  the  perpendicular 
SA,  the  point  A  becomes  a  source  of  disturbance  which  is  pro¬ 
pagated  after  this  moment  into  both  media.  According  to  Huygens 
principle  each  disturbed  element  or  point  is  to  be  considered  as 
a  source  generating  the  elementary  waves.  For  example,  the 
reflected  wave  ^  is  the  envelope  of  the  system  of  halfspheres  having 
centres  inside  the  circle  of  radius  AB  in  the  plane  z=.Q  (Fig.  1). 

If  the  velocity  of  propagation  in  the  second  medium  (  v^  )  is  larger 
than  that  in  the  first  medium  (  v^  ),  the  corresponding  half  spheres 
in  the  second  medium  will  have  larger  radii  and  the  front  of  the 
refracted  wave  will  be  more  advanced  than  that  of  the  reflected. 
Therefore,  the  points  B',  B"  and  C  of  the  ring  BC  (  and  certain 
points  nearer  to  A  )  will  become  new  centres  of  disturbance  before 
or  at  the  moment  the  direct  wave  has  reached  the  circle  B.  Thus, 
the  other  front  wave,  say  TGf,  should  be  formed  in  the  first  medium 
by  disturbances  generated  at  the  points  of  a  ring,  say  KC. 
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To  simplify  the  discussion  we  assume  now  that  there 

are  only  compre s sional  waves.  The  disturbance  starts  to 

propagate  in  the  second  medium  at  the  moment  t^^h/v^  , 

at  which  the  incident  wave  reaches  the  point  B'  (Fig.  3). 

A 

Let  ~t  be  the  time  necessary  for  the  incident  wave  to  travel  over 

a 

the  path  B'B;  £  B'B/v^  .  During  the  same  time  the  dis  - 
turbance  will  propagate  through  the  second  medium  in  all 

A 

directions  and  will  reach  a  point  K  where  the  distance 
Whether  AK=  B 'B  depends  on  the  velocity  .  However,  if 
v2  ^  V1  *  ^  'B  <  AK.  Since  B  'B  /AB  tends  to  zero  if  B '  tends  to 
A,  we  can  always  find  a  point  B  which  will  be  reached  by  the 
incident  wave  earlier  than  by  the  refracted.  On  the  other  hand 
there  are  obviously  points  P  for  which  the  time  of  arrival  of  the 
direct  wave  is  larger  than  that  of  the  refracted  or,  more  exactly, 
of  the  disturbance  propagated  in  the  second  medium.  Thus,  at 
each  moment  the  interface  may  be  divided  into  three  regions: 
the  central,  in  which  the  first  disturbance  comes  from  the  first 
medium,  a  ring  in  which  the  points  become  centres  of  disturbance 
through  the  second  medium  in  first  and  the  outer  region  which  is 
at  rest. 

We  can  now  find  the  boundary  between  the  first  two  regions 
represented  by  the  point  M  (Fig.  3)  by  the  condition  that  the  direct 


wave  from  M'  and  the  refracted  from  C  will  reach  M  at  the  same 


moment.  Putting  and  ^^s^our  condition  is 

MM'  /CM  —  ^  .  From  the  triangle  CMM'  it  follows  that 


MM'sCM  sin^C.  MCM'  and  ,  on  neglecting  a  small  difference, 


we  can  write  M' CM=  90 -  O'  ,  where  $  is  the  angle  of  incidence. 

,  is  the  angle  of  total  reflection  and  the 
circumference  of  the  radius  fl  £c*n  &is  the  boundary  referred  to 


earli  er . 

To  obtain  the  wave  front  of  refractions  we  now  have  to 
take  into  account  the  fact  that  the  points  of  the  interface  become 
sources  at  different  times,  since  the  disturbance  is  propagated 
as  concentric  circles  from  the  point  A.  In  order  to  write  the 
equation  of  spherical  waves  emitted  by  one  of  these  points,  say 
C,  we  denoted  by  r  ,  the  polar  coordinates  of  the  center  of 
these  waves.  The  radius  of  the  sphere  is,  growing  in  the  upper 
medium  with  a  velocity  lFf  and  is  equal  to  l£^VZ^)if  ^denotes  the 
time  when  C  (Fig.  4)  starts  to  emit  a  disturbance.  If  we  put  — 
Zc  -f-  and  T-~t~Co the  equation  of  the  spherical  wave  is 


To  find  the  envelope  of  these  surfaces  we  have  to  form 

the  eqi 


or 


(4) 
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Hence 


■Un.  Cfc~l  , 


and  by  (5)  we  obtain 


q 7  =  7==-  )  Q>c 

'c  SZ+Y'-  7 


y 


=-  (6) 


\scx+  yz 


_  r. 


cc -zcc*<fc  =  if(r-v 

2.  2. 


y-zt**<fc=  z£(r-  Yc 

Now,  on  inserting  (6)  in  (5)  it  is 

f  /  -  £ )  7«  *  ~  % rr 


(7) 


(8) 


and  by  (7)  and  (3) 

V 


2 


or 


(V- 1*}  si  (r- 1  f+  z' 


X  x  ^ 

cj 

i 

Since  we  consider  the  case  it  is 

I  *? 


•  §•  fr-  ^ 1  z 


(9) 


On  eliminating  /-from  (8)  and  (9)  we  obtain 


2=  - 


^/»i 


t =====-- V  ^  ^ 1 -  ^=i^7<10) 

\T  -  ^  V  ?  x/7^r>  *  * 

/  ,Z 

It  is  Vi fy^.  $M $ %  Therefore,  on  taking  the  sign  -f  if  z  is 


positive  in  the  upper  medium,  we  have 


(11) 


This  equation  represents  the  envelope  of  elementary  Huygens 
waves  which  has  the  form  of  a  cone.  The  generatrix  of  this  cone  is 
inclined  to  the  interface  at  the  angle  of  total  reflection  0  .  If  we 


consider  the  case  of  a  source  located  at  the  interface  (Fig  4)  and  the 
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radius  of  the  direct  wave  in  the  upper  medium  at  an  instant 
is  Vft ,  it  may  be  easily  seen  from  this  figure  and  equation  (11) 
that  the  envelope  is  tangent  to  the  direct  wave  W^.  It  is  also 
obvious  that,  if  the  source  is  above  the  interface,  the  reflected 
wave  will  play  the  part  of  the  direct  in  the  former  case  and  the 
envelope  will  be  tangent  to  it. 

By  this  elementary  representation  the  physical  conditions 
of  formation  of  refractions  as  well  as  the  form  of  the  wave  and  its 
propagation  are  clarified  and  it  is  possible  to  draw  some  further 
conclusions. 

We  have  seen  that  the  waves  of  the  'fourth"  kind  exist  in 
the  upper  medium  under  the  assumption  that  the  latter  is  of  lower 
velocity.  These  waves  are  generated  successively  by  points  of  the 
interface  because  of  the  faster  propagation  of  the  disturbance  in  the 
second  medium.  Moreover,  these  waves  start  at  the  points  of  a 
circle  at  which  the  incident  wave  in  the  first  medium  impinges  at 
an  angle  equal  to  the  angle  of  total  reflection,  i.  e.  for  a  certain 
r  >0^  and  being  determined  for  larger  values  of  r,  have  the  simple 
form  of  a  regular  frustrum.  In  order  to  avoid  any  conf  usion  a  name 
should  be  given  to  these  waves  which  would  correspond  to  their 
nature.  Perhaps  "interface  generated  waves"  is  the  most  appropri¬ 
ate  name  for  these  waves  whose  apparent  velocity  parallel  to  the 
interface  is  V2  and  whose  velocity  normal  to  the  wave  front  is  v^  , 
the  wave  front  being  inclined  at  angle  0  from  the  interface. 
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We  have  not  yet  taken  into  account  the  ability  of  a  medium 
to  transmit  longitudinal  and  transverse  waves.  To  generalize  our 
former  result  we  have  to  assume  that  if  a  disturbance  is  transmitted 
in  a  second  elastic  medium  the  points  of  the  interface  become  sources 
of  both  kinds  of  waves.  Simple  considerations  lead  u&  then  to  the 
conclusion  that  there  must  exist  those  kinds  of  waves  which  are 
represented  by  the  Figure  2. 

It  is  also  important  that,  according  to  the  above  representation, 
the  number  of  "rays"  concurrent  at  a  certain  point  of  the  interface 
vary  with  the  position  of  the  point.  Let  us  consider  the  simplest 
case  where  there  are  only  compression^  waves  (Fig.  1).  If  we 
imagine  one  pulse  propagating  as  a  spherical  wave  from  the  source 
above  the  interface,  there  will  be  three  rays  concurrent  simultaneously 
at  each  point  inside  the  circle  first  reached  by  the  direct  wave,  since 
the  disturbances  from  other  points  of  the  interface  travelling  through 
the  second  medium  come  too  late.  At  the  point  M  (Fig.  3,  corres¬ 
ponding  to  the  angle  of  total  reflection)  this  disturbance  comes  at  the 
same  moment  as  the  direct  wave  and  the  Schmidt's  "fourth"  ray  will 
coincide  with  the  reflected  ray.  Thus,  inside  of  this  region  there 
are  concurrent  three  rays,  as  it  was  always  assumed,  but  in  the 
outer  region,  where  the  refracted  wave  reaches  the  interface  in  the 
second  medium  earlier  than  the  direct,  only  interface  generated 
wave  is  present  in  the  first  medium  and  only  two  "rays"  are,  there- 
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fore,  concurrent  at  each  point.  The  Schmidt's  remark  that 
there  may  be  at  one  point  of  interface  four  rays  concurrent 
(due  to  the  same  pulse,  which  fact  would  contradict  the  laws  os 
Optics)  does  not  hold. 

The  above  representation  of  interface  generated  waves  may 

be  generalized  for  the  case  of  several  layers.  There  should  be 

n 

waves  generated  at  all  interfaces.  The  wave  ,  for  example, 
generated  at  the  second  interface  will  be  refracted  at  the  first 
interface  and  enter  the  first  medium.  A  receiver  at  the  free 
surface  might  register  this  wave  as  a  first  arrival.  This  will 
occur  if  l^<l^<Z^and  if  it  is  placed  at  a  distance  far  enough  from 
the  source  S.  The  region  which  was  in  case  of  two  media  subjected 
to  the  action  of  a  single  wave  coming  through  the  second  medium  will 
be  crossed  by  waves  coming  from  the  third  medium  and  the  whole 
picture  becomes  more  complicated.  But  as  to  the  first  arrivals, 
it  is  also  obvious  that  there  will  be  a  limit  distance  from  the  source 
at  which  the  first  arrival  is  the  direct  wave  and  that  at  a  larger 
distance  the  first  arrival  corresponds  to  some  interface  generated 
wave  if  the  velocity  of  propagation  in  the  first  medium  is  not  the 


largest  one. 
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CONCLUSIONS 


The  displacements  representing  a  disturbance  propagating 
in  one  or  several  elastic  media  are  usually  expressed  in  terms 
of  solutions  of  wave  equations.  If  we  write  a  general  solution  of 
a  wave  equation  in  the  form  which  follows  from  that  given  by 
Whittaker,  all  known  solutions  corresponding  to  particular  cases 
may  be  derived  from  it.  Thus,  it  was  shown  that,  if  a  most 
general  linear  form  of  the  variables  x,  y,  z  and  t  is  taken  as  an 
argument  of  an  arbitrary  function  by  means  of  which  the  general 
solution  is  built  up,  the  so-called  steady  state  solution  as  well 
as  those  used  by  Coulomb,  Schermann,  Cagniard  and  others  are 
obtained  by  specification  of  conditions. 

In  Part  II  the  generation  of  the  so-called  '"refraction"  waves 
is  represented  in  an  elementary  way.  This  theory  clarifies  the 
physical  conditions  leading  to  the  "fourth"  kind  of  waves,  the 
direct,  reflected  and  refracted  ones  being  the  first  three.  It  shows 
also  the  different  boundary  conditions  for  different  regions  of  an 
interface. 

In  concluding  this  paper  the  author  wishes  to  thank  Professors 
Maurice  Ewing  and  Frank  Press  for  the  very  valuable  discussion  of 


this  work. 
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Fig.  I  Formation  of  wave  fronts 

(Source  above  the  interface) 


Fig.  2  A  complete  set  of  waves  in  two  solid  media 
produced  by  a  point  source,  according  to  Cagniard  (1939). 
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Fig.  3  Different  regions  reached  by  waves 


Fig.  4  Formation  of  wave  fronts 
(Source  at  the  interface) 
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